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ABSTRACT 

A candidate accreting massive black hole (BH) with Mbh ^ 10^ Mq has recently been identified 
at the center of the dwarf starburst galaxy Henize 2-10 (He 2-10). This discovery offers the first 
possibility of studying a growing BH in a nearby galaxy resembling those in the earlier universe, and 
opens up a new class of host galaxies to search for the smallest supermassive BHs. Here we present 
very long baseline interferometry observations of He 2-10 taken with the Long Baseline Array (LBA) 
at 1.4 GHz with an angular resolution of ~ O'.'l x 0^.^03. A single compact radio source is detected at 
the precise location of the putative low-luminosity active galactic nucleus. The physical size of the 
nuclear radio emission is < 3 pc x 1 pc, an order of magnitude smaller than previous constraints 
from the Very Large Array (VLA), and the brightness temperature of Tb > 3 x 10^ K confirms a 
non-thermal origin. These LBA observations indicate that the nuclear radio emission originates from 
a single object and exclude the possibility of multiple supernova remnants as the origin of the nuclear 
radio emission previously detected with the VLA at lower resolution. A weaker, more extended, off- 
nuclear source is also detected with the LBA and a comparison with mult i- wavelength ancillary data 
indicate that, unlike the nuclear source, the off-nuclear source is co-spatial with a super star cluster, 
lacks a detectable X-ray point-source counterpart, and is almost certainly due to a supernova remnant 
in the host star cluster. 

Subject headings: galaxies: active — galaxies: dwarf — galaxies: individual(He 2-10) — galaxies: 
starburst 



1. INTRODUCTION 

Supermassive black holes (BHs) with Mbh ^ 10^ — 
lO^M© are thought to reside in t he nuclei of es sentially 
all m assive galaxies with bu lges (Kormendy fc Richstone| 
1995 Magorrian et al.|1998 ), however the origin of these 
BHs is largely unknown. Dwarf galaxies with low masses 
and relatively quiet merger histories are potential hosts 
of the least-massive BHs, and can therefore provide valu- 
able constraints on the properties of the first primordial 
^^seed" BHs (e.g. [ Volonteri et al. 2008| [van Wassenhove 

~Bl 



et al.|2010[[Bellovary et al |201l[ )." The IJH mass distribu- 
tion and occupation fraction in modern dwarfs can help 
distinguish between various seed formation mechanisms 
at high-redsh ift, such as stellar niergers in compact star 



clusters (e.g. jGiirkan et al. [2004 Devecchi 

)puL 



20091), remnants from Population III stars (e.g. MacLau 



Volonteri 



& Rees||2001l), or the collap se of mass ive objects from 



proto -galactic gas (e.g. [ Bromm fc Loeb || 2Q03[ |Begelman| 
et al.j,20Q6 . ^Lodato fc Nat arajan| 2006 [ ) Active galactic 
nuclei (AGN) in dwarf galaxies undergoing significant 
star formation are particularly important for studying 
BH growth in physical conditions similar to those ex- 
pected in the earlier universe. 

Observationally, however, few dwarf galaxies are 
known to host massive (as opposed to stellar-mass) BHs. 
The most well-studied systems in the low-mass regime 
are the low-luminosity active galactic nuclei (LLAGN) 
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in NGC 4395 ( [Filippenko fc Sargent [[19891 [Fihppenko fc 
Ho|2003 ), an essentially bulgeless late- type spiral galaxy 
witn a nuclear star cluste r, and Pox 52 (Barth et al. 
Thornton et al.|20'Q8 ), a dwarf elliptical. Estimates 
le BH masses in th ese galaxies a re ^ few x IO^M^t) 
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(iPilippenko fc Ho[[2003 



Barth et al. 



al. 20 65|. Systematic searche s for ot 
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Peterson et 
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by Greene fc Ho[ ( [2004 , [2007[ ) have revealed - 200 BHs 
with virial masses of Mbh 10^ — lO^M© in a sam- 
ple of broad-line AGN from the Sloan Digital Sky Sur- 
vey. Although most of the host galaxies are larger and 
brighter than typical dwarfs (the vast majority have ex- 
tended disks), they are sub- L* galaxies witho ut classical 
bulges (jGreene et al. 2008; Jiang et al.||2011|). 



Recently, Reines et"ar[ ^2011) have discovered a can- 
didate accreting massive BH at the center of Henize 2- 
10 (He 2-10), a nearby (9 Mpc), gas-rich, dwarf galaxy 
in the mid st of an intense burst of star formation (e.g. 
Allen et al , 1976; Johansson 1987] [Kobulnicky et al.|1995l 



Johnson fc Kobulnicky 2003i). He 2-10 has a compact 
1 kpc) irregular morphology (Figure [T]) with no dis- 
cernible bulge, and hosts a large population of young, 
massive and dense "super star clusters" that are thought 
to be progenitors of globular clusters dConti fc Vacca 



1994[ [Johnson et al. 2000; Chandar et al. 2003). In fact, 
the study of i nfant massive star clusters (Reines et al. 
2008a|b[ 120101 ) in He 2-10 is what ultimatel y led to the 
serendipitous discovery of the putative BH ( [Reines et al. 



201 1[ ). 

The LLAGN is detected as a central unresolved radio 
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Fig. 1.— T hree-color HST imag e of He 2-10 with 8.5 GHz VLA 
contours from lReines et ar] ( |201l{ . The dashed box indicates the 
field of view or the LBA image shown in Figure p] 



source with the Very Large Array (VLA^; Kobulnicky 



al. 



jQhnsQn|[T999l IJohnson Kobulnicky||2UU3| [Keineset 



2011) that is spatiahy coincident with a hard X-ray 



point source observed with the Chandra X-ray Observa- 



tory (lOtt et al ||2QQ5| [Kobulnicky &: Martin|2QlQ 
et ai.[201l[ ) and clearly separated from the promii 



Reines 

et al.|201ip and clearly separated from the prominent Hn 
regions m the galaxy. No star cluster is seen at the loca- 
tion of the nuclear radio source even in deep near-infrared 



Hubb le Space Telescope (HST) observations ( [Reines et al 



2011), and the position of the radio source is consistent 
with the dynamical center of the galaxy. Moreover, the 
nuclear radio source lies at the center of a ~ 250 pc-long 
ionized gas structure with a velocity gradient that is con- 
sistent with an outflow or rotating disk of ionized mate- 
rial. The relative strength (i.e. the ratio) of the radio 
to hard X-ray luminosities of the c ompact ce ntral source 
is consistent with known LLAGN ('Ho""2008), but differs 
by orders of magnitude from active X-ray binaries and 
supernova remnants (SNRs). A young supernova origin 
is ruled out since the VLA flux density at 8.5 GHz has 
not measura bly declined in more than a decade (from 
1994 to 2004; [Johnson fc Kobulnicky|[2QQ3l [Reines et aC 
2011) and Very Long Baseline Array (VLBA) observa- 
tions at 5 GHz do not detect any sources on scales of 
0.1 pc (a typical size for a young compact radio super- 

The VLBA non-detection 
however, does not rule out 
Using an empirical rela- 



nova; 
of the 



Ul vest ad 
nuclear 



et al. 



2007) 



source, 

a more extended LLAGN jet. 
tionship between compact radio emission, compact hard 
X-ray emission and BH mass (the ^^fundam ental plane of 



black hole activity," Merloni et al. 2003), Reines et al. 
(2011) estimate a mass of \og[MBYi/ Mq) = 6.3 ± 1.1 for 
the BH in He 2-10. 

As part of an effort to verify the nuclear source in 
He 2-10 is indeed a LLAGN, we are undertaking follow- 
up observations of the nuclear radio emission using very 
long baseline interferometry (VLBI). In this Letter, we 
present VLBI observations of He 2-10 taken with the 
Long Baseline Array (LBA) at 1.4 GHz that have a beam 
solid angle ^ 38 x s maller than the VLA observations of 
Reines et aL]([2011[) and - 128 x larger than the VLBA 
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Fig. 2. Top: Final cleaned LBA image of He 2-10. The brighter 
central source is the candidate LLAGN. The beam (106 mas x 
29 mas) is shown in the lower left corner. Bottom: A close-up view 
of the two detected sources. Contour levels are 4,8 and 16 x the 
rms noise (Icr = 32/iJy/beam). The colorbar has units of Jy/beam. 



observations of 'Ulvestad et al. (2007) that apparently 
resolved out the nuclear source. Our LBA observations 
reveal a parsec-scale radio source coincident with the can- 
didate LLAGN, supporting the case for a massive BH in 
He 2-10. 

2. OBSERVATIONS AND DATA REDUCTION 

We observe d He 2-10 on 2011 Ju ly 22 with the LBA in 
eVLBI mode ( [Phillips et al.|[2QQ7 ). Continuum observa- 
tions at 1.4 GHz (21 cm) were made using the four eVLBI 
capable antennas at this frequency - the Australia Tele- 
scope National Facility (ATNF) observatories at Parkes 
(64-m single dish) , the Australia Telescope Compact Ar- 
ray (ATCA; five 22-m antennas phased), and Mopra (22- 
m single dish), as well as the Hobart antenna (26-m single 
dish) operated by the University of Tasmania. The ob- 
servation duration was 8 hours, but the target source was 
only visible to Parkes for 5 hours, due to the relatively 
restrictive Parkes elevation limit. Two bands of width 
64 MHz were sampled with two bit precision in the fre- 
quency ranges 1.368 - 1.432 GHz and 1.432 - 1.496 GHz 
at the three ATNF stations, while only the lower fre- 
quency band was sampled at Hobart due to the more 
limited network connectivity to this station. The data 
(with an aggregate rate of 3.5 gigabits per second) were 
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TABLE 1 
Single Gaussian Model Fits 



Parameter 


Nuclear source 


Off-nuclear source 


Right Ascension (J2000) 


08^36™15.1200^ ± 0.0003 


08^36™15.2016^ ± 0.0004 


Declination (J2000) 


-26°24'34.15r' ± 0.004 


-26°24'33.732'' ± 0.005 


Integrated flux density (mJy) 


0.98 ± 0.06 ± 0.20 


0.59 ± 0.08 ± 0.12 


Peak flux density (mJy/beam) 


0.73 ± 0.03 ± 0.15 


0.29 ± 0.03 ± 0.06 


Fitted Major axis (miUiarcseconds) 


124 ± 5 


137 ± 13 


Fitted Minor axis (miUiarcseconds) 


33 ± 1 


44 ± 4 


Fitted Position angle (degrees) 


77t\ 


71+3 


Deconvolved Major axis (miUiarcseconds) 


65 ± 10 


90 ±20 


Deconvolved Minor axis (miUiarcseconds) 


10 ± 10 


20 ±20 


Deconvolved Position angle (degrees) 


68ll 





Note. — Flux densities are quoted as fitted value ib fit error it absolute flux density scale error (estimated at 20%). At 
the adopted distance to He 2-10 (9 Mpc), 10 mas 0.44 pc. 



correlated on a computer cluster at the ATCA in real 
time using the DiFX software correlator (Deller et al. 
2QQ7D 



The observations utilized a phase reference cadence 
of 5-minute loops with 3.5 minutes on-source and 1.5 
minutes on a phase calibrator separated by 2.4° (J0846- 
2607). Additionally, a fringe finder (0537-441 or 0834- 
291) was observed for 5 minutes before and after every 
^2-hour block on He 2-10 and the phase calibrator. Af- 
ter editing, a total of ~3.5 hours of usable integration 
time was obtained on He 2-10, of which the first 50 min- 
utes lacked the Parkes antenna. The predicted sensi- 
tivity for the target field using natural weighting was 
28 /iJy/beam. However, data from the phase reference 
source showed that the phased ATCA station suffered 
from severely degraded sensitivity during the first half 
of the experiment (reduced by a factor of ~3). Instru- 
mental failure or poor phasing of the array are potential 
explanations, but the actual cause could not be deter- 
mined based on system logs and monitor data. As the 
ATCA is the second most sensitive station in the array, 
this corresponds to a significant (10%) reduction in the 
predicted final image sensitivity to 31 /iJy/beam. 

Data reduction was performed using the AIPS^ soft- 
ware package, an d imaging was performed using Difmap 
([Shepherd' 1997). The data were first amplitude cali- 
brated using the logged system temperature data, de- 
lay calibrated using the phase reference source, and then 
phase and amplitude calibration were performed using 
a model of J0846-2607 derived from these observations. 
The absolute amplitude calibration accuracy of the visi- 
bility data is expected to be at the 20% level. The target 
field data were split, averaged in frequency and exported 
from AIPS in uvfits format for imaging in Difmap. 

In Difmap, a 1024 x 1024 pixel image with a pixel size of 
5 miUiarcseconds (mas) was formed using natural weight- 
ing. Two sources were detected and cleaned. After clean- 
ing, a phase self-calibration was applied with a solution 
interval of 15 minutes, the shortest time at which reliable 
solutions could be obtained. This self-calibration step in- 
creased the recovered flux density of the two sources by 
~40%. Figure [2] shows the final cleaned image. The la 
image sensitivity is 32 /iJy/beam, consistent with expec- 
tations, and the angular resolution is 106 mas x 29 mas, 

^ http:/ /www. aips.nrao.edu/index.shtml 



at a position angle of 80°. At the adopted distance of 
He 2-10 (9 Mpc), this angular resolution corresponds to 
a linear size of ~ 4.6 pc x 1.3 pc. 

To facilitate the comparison of multi-wavelength im- 
ages of He 2-10 , we registe r our LBA image and the data 
presented in Reines et al. (2011) to the same astromet- 
ric reference frame. We make the reasonable assumption 
that the two LBA sources coincide with two compact, 
non-thermal VLA sources to derive the astrometric shift. 



The sources injReines et al. (2011) have positions refer- 
enced to the Two Micron All Sky Survey catalog ( jSkrut-] 
skie et al.|2006D and are offset by - 0'.'04 west and - 0709 
north relative to the LBA positions. 

3. RESULTS 

We fit each of the two detected sources with a single 
Gaussian model using the task JMFIT in AIPS, the re- 
sults of which are shown in Table [H 

3.1. The Nuclear Source: LLAGN 



The candidate LLAGN identified by [Reines et al.| 
(2011) is clearly detected in our LBA observations at 
1.4 GHz and is compact on parsec scales. Although the 
source appears marginally resolved, the significance of 
this measurement is low. In addition, residual phase and 
amplitude calibration errors will certainly remain after 
phase referencing and self-calibration, since the atmo- 
spheric phase will vary on timescales shorter than the 
solution interval of 15 minutes, and the signal-to-noise of 
the calibration solutions is relatively low. These residual 
errors will have the effect of broadening the point-source 
response, and so we consider the deconvolved source size 
(plus errors) to be an upper limit on the actual source 
size. This upper limit of 75 mas x 20 mas corresponds 
to a linear size < 3 pc x 1 pc. C ompared to the VLA 



observations of Reines et al. (2011), which have a linear 
resolution of ~ 24 pc x 9 pc, this yields an order-of- 
magnitude reduction in the maximum source size. 

The lower-resolution VLA observations allowed for the 
possibility of multiple discrete sources contributing to 
the flux densit y in the nuclear region, such as a collec- 
tion of SNRs ( juivestad et~aL]|2007[ [Reines et an|2011 ) 
with each remnant having a t ypical size of ^ tew pc 



(e.g. M82; iFenech et al.[ [2010[ ). However, we only de- 
tect one parsec-scale source at the center of this region 
(at the ~ 31<j level), despite the LBA observations hav- 
ing a beam solid angle ~ 38x smaller than the VLA 
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Fig. 3. — Left: HST Pao; image from 'Reines et al.'("2011') of the central region of He 2-10 showing ionized gas emission. Radio contours 
are in green (LBA) and black (VLA) and the beam sizes are shown in the lower right corner. The candidate LLAGN is at the center of 
the ~ 250 pc-long ionized gas structure. Right: HST F814W (/-band) image of the same region showing young super star clusters. The 
off-nuclear LBA source is spatially coincident with the most luminous stellar cluster and is almost certainly related to a SNR. 



observations. This indicates the nuclear radio emission 
comes from a single object and argues against multiple 
SNRs within the unresolved nuclear VLA source. It is 
also very unlikely that multiple SNRs reside within the 
nuclear LBA parsec-size source, given that VLBI obser- 
vations of M82 reveal dozens of well-resolved SNRs that 
are clearly separated from one an other and spread o ver 



the central kpc of the galaxy (e.g. Fenech et al. 

Using the fitted peak flux density of 0.73 m^ 
and the synthesized beam size of 106 mas x 29 mas, we 
calculate a lower limit to the brightness temperature of 
Tb > 1.4 X 10^ K. However, as described above, the emit- 
ting region is likely signifi cantly smaller than the beam 
size. Therefore, following Blundell et al. ( 1996 ), the lower 
limit on the brightness temperature can be increased to 
Tb > 2.9 X 10^ K using the upper limit on the source size 
of 75 mas x 20 mas. This high brightness temperature 
rules out free-free emission from thermal electrons with 
a characteristic temperature of ~ 10^ K (i.e. an Hll re- 
gion), and confirms a non-thermal origin for the nuclear 
radio source as indicated by the spectral index derived 
from VLA observa tions at higher frequencies and lower 
angular res olution ( Johnson fc Kobulnicky|[2QQ3 Reines 
eFaI^[2QTT] ). 



The mtegrated flux density of the nuclear LBA source 
is 0.98 dz 0.21 mJy. At a distance of 9 Mpc, this cor- 
responds to a spectral luminosity of Li.4ghz ^ 10^^ W 
Hz-^ For comparison, this is approximately an order 
of magnitude more luminous than the LLAGN in NGC 
4395 at similar spatial scales (i.e. intermediate of VLA 
and VLBA obs ervations of NGC 4395; Figure 2 in|^^ 
|bel et al.|[2QQT] ). 

I'he measured flux density is roughly consistent with 
predictions based on higher frequency VLA observations 
^Reines et al. 2011). We have re-derived the VLA flux 
densities using the AIPS task JMFIT to be consistent 
with the method used here, fitting a point source plus 
a background (level + slope) to the region immediately 
surrounding the nuclear source. We obtain flux densities 
of 0.86±0.02 and 0.64±0.01 mJy at 4.9 and 8.5 GHz, re- 
spectively. Tl]^s£_vBhiesax^ lower than those 
quoted in Reines et al. (2011), which were derived using 
larger apertures and ditterent background regions (see 
the Supplementary Information). We note, however, that 



this does not affect the results of that paper in any ap- 
preciable way. Assuming an additional 15% uncertainty 
for the absolute flux calibration of each VLA image, the 
spectral index a = —0.54 ± 0.39 {Sj^ ex u'^). The pre- 
dicted 1.4 GHz flux density is then 1.7^o ^ mJy. Thus, 
the integrated LBA 1.4 GHz flux density of 0.98 ± 0.21 
mJy is (marginally) consistent with expectations if the 
spectral index remains constant between 1.4 and 8.5 
GHz. However, our data do not exclude the possibility 
that additional flux is present at 1.4 GHz on spatial scales 
intermediate between the LBA and VLA beams, similar 
to what is observed in steep-spectrum Seyfert nuclei in 
which the radio emission is thought to be domi nated by 
low-surface brightness features such as jets |Orienti fc 
^Prieto^^2010J . It is also possible that the radio spectral 
energy distribution flattens at lower frequencies and ex- 
hibits a spectra l turnover, a common c haracteristic of 
AGN cores (e.g. [Pearson et al.|p92| [OT) ea 1998). 



3.2. The Off- Nuclear Source: Supernova Remnant 

A second source is detected in our LBA observations 
that is located to the north-east of the nucleus (Fig- 
ure [2|. This source is fainter and more extended than 
the nuclear source, with Tb ~ 4.0 x 10"^ K. The bright- 
ness temperature of the nuclear source is therefore ~ 7x 
larger than the off- nuclear source. Additionally, in con- 
trast to the nuclear source, the off- nuclear s ource lacks 
a detectable X-ray point-source counterpart (lOtt et al 



2005[|Kobulnicky fc Martin|2010| [Reines et al.|201ip and 

IS spatially coincident with a lummous super star clus- 
ter (Figurejsl), strongly suggesting the radio emission is 
from a SNK in the host cluster. The most massive O 
stars in a super star cluster are expected to begin ex- 
ploding as supernovae at an age of ~ 3.5 Myr (Leitherer 
jCt al.|1999| ) and the host star cluster is the most luminous 
in the gal axy with a mass of ^ 2 x 10^ Mq and an age of 
~ 5 Myr fChandar et al. 2003). For comparison, the 1.4 
GHz radio luminosity of the off- nuclear LBA source in 
He 2-10 is about twice that of the Galactic SNR Gas A, 
assuming the flux den sity and spectral index measured 
by Baars et al. (1977) and the distance derived by Reed 
|et dA.\(\ld%l 
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4. CONCLUSIONS AND DISCUSSION 

We observed He 2-10 with the LB A at 1.4 GHz with 
an angular resolution of ~ 0^1 x 0V03, corresponding to 
a linear size of ~ 4.6 pc x 1.3 pc. Our findings are 
summarized below: 

1. We clearly detect compact radio e mission from 
the candidate LLAGN identified by Reines et al. 
(^2011), supporting the case for an accreting mas- 
sive BH in this dwarf starburst galaxy. The physi- 
cal extent of the nuclear radio emission is <3 pc x 
1 pc, an order of magnitude smaller than previous 
constraints from the VLA. 

2. Our observations indicate the nuclear radio emis- 
sion originates from a single object. The possibility 
of multiple SNRs as the origin of the nuclear radio 
emission previously detected by the VLA on larger 
scales is disfavored by the detection of only one 
bright compact LBA source in the nuclear region. 

3. The high brightness temperature of the nuclear 
source, > 3 x 10^ K, rules out an Hll region 
and confirms a non-thermal origin. 

4. The integrated fiux density of the nuclear LBA 
source at 1.4 GHz is ~ 1 mJy, corresponding to 
a spectral luminosity of Li.4ghz ^ 10^^ W Hz~^. 
This is an order of magnitude larger than the 
LLAGN in NGC 4395 at similar spatial scales. 

5. We detect a weaker and more extended off- nuclear 
source that, unlike the nuclear source, lacks a 
detectable X-ray point-source counterpart, is co- 
spatial with a super star cluster, and is almost cer- 
tainly due to the emission from a SNR in the host 
star cluster. 

While multi-wavelength observations provide a com- 
pelling c ase for the presen ce of a LLAGN in He 2-10 (see 
^andjR emes et al. 2011L even radio observations alone 



2007D and the 
.5 GHz observations from the 



are placing stringent constraints on the nuclear source. 
An Hll region is ruled out by the high brightness temper- 
ature derived from the LBA observations presented here 
and the spectr al index derived from previous VLA obser - 
vations ( [Johnson fc Kobulnicky|2003 Reines et al.|2011 ). 
A youn^ compact supern ova is ruled out by the VLBA 
non-detection at 5 GHz (jUlvestad et al. 
lack of fiux evolution in 
VLA. Multiple SNRs are ruled out by our LBA obser- 
vations that detect only one parsec-scale nuclear source. 
This leaves only two remaining possibilities for the nu- 
clear radio emission: a LLAGN or a single SNR at the 
center of the galaxy, far from any detectable signs of star 
formation. The luminosity of the nuclear source would 
be high for a SNR - it is more powerful than 31 out of 
36 known radio S NRs in M82 - and it would likely be 
just decades old (iFenech et 
tions with the LBI 



2010| ). Future observa- 



S Fenech et al.| 
I (at higher frequency) and the VLBA 
(at low frequency) can be employed to further charac- 
terize the nuclear radio source in He 2-10, providing a 
high-resolution radio spectral energy distribution and de- 
tailed morphological information necessary for a secure 
identification. 
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